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ABSTRACT 

We have obtained multi-object spectroscopy of H II regions in the spiral galaxy M5 1 with the Keck I telescope 
and the Low Resolution Imaging Spectrometer. For ten objects we have detected the auroral line [N II] A5755, 
while [S III] A6312 has been measured in seven of these. This has allowed us to measure the electron temper- 
ature of the gas and to derive oxygen, sulfur and nitrogen abundances for the ten H II regions. Contrary to the 
expectations from previous photoionization models of a few H II regions in M5 1 and from strong-line abun- 
dance indicators, the O/H abundance is below the solar value for most objects, with the most metal-rich H II 
regions, P203 and CCM 72, having log (O/H) = -3.16 [~1.4(O/H) ] and log (O/H) = -3.29 [~1.O(O/H) ], 
respectively. The reduction of O/H by factors up to two or three with respect to previous indirect determinations 
has important consequences for the calibration of empirical abundance indicators, such as R-a, in the abundance 
and excitation range found in the central regions of spiral galaxies. The abundance gradients in these galaxies 
can therefore be considerably flatter than those determined by making use of such empirical calibrations. The 
Hll regions with a measured electron temperature span the 0.19-1.04 Rq range in galactocentric radius, and 
indicate a shallow abundance gradient for M51: -0.02 ±0.01 dexkpc -1 . The S/O abundance ratio is found to 
be similar to previous determinations of its value in other spiral galaxies, log (S/O) w— 1.6. Therefore, we find 
no evidence for a variation in massive star initial mass function or nucleosynthesis at high oxygen abundance. 
An overabundance of nitrogen is measured, with log (N/O) ~ -0.6. Based on our new abundances we revise 
the effective yield for M51, now found to be almost four times lower than previous estimates, and we discuss 
this result within the context of chemical evolution in galactic disks. Features from Wolf-Rayet stars (the blue 
bump at 4660 A and the C III line at 5696 A) are detected in a large number of H II regions in M51, with the 
C III A5696 line found preferentially in the central, most metal-rich objects. 

Subject headings: galaxies: abundances — galaxies: ISM — galaxies: spiral — galaxies: individual (NGC 
5194) — H II regions 



X 



1. INTRODUCTION 

Measuring the chemical abundances of gaseous nebu- 
lae is a crucial step to understand the chemical evolu- 
tion and the nucleosynthesis in spiral galaxies. Histori- 
cally, one of the main obstacles has been the unavailabil- 
ity of electron temperatures (T e ) at high metallicity. This 
is due to the enhanced cooling via far-IR lines causing 
the faint auroral lines that are necessary to determine T e 
(e.g. [Om]A4363, [Nn]A5755, [S ni] A6312) to drop be- 
low detectability in the spectra of Hll regions. This dif- 
ficulty affects the observations of nebulae lying in the in- 
ner portions of most spiral galaxies, where the oxygen 
abund ance can exceed the solar va l ue flo g (O/H) =-3.31, 
lAllende Prieto. Lambert & Asplundl 1200 lH . In order to 
make this important astrophysical problem tractable, strong- 
line methods have been calibrated and almost univer- 
sally adopted for the study of abundance gradients in spi- 
rals and, more recently, to estimate ch emical abundances 
in high-redshift star-forming galaxies dAlloin et.aU [l979; 

1 The data presented herein were obtained at the W.M. Keck Observatory, 
which is operated as a scientific partnership among the California Institute of 
Technology, the University of California and the National Aeronautics and 
Space Administration. The Observatory was made possible by the generous 
financial support of the W.M. Keck Foundation 



iPagel et alJ Il979t Edmunds & Pagel 1984; Dopita & 
Evans 1986: lMcGaughlll991t Diaz & Perez-Montero 2000; 
| Pilvuginl2000HKewlev & Dopital2002t iDenicolo et al.ll2002t 
Kobulnick v. Henrv & Phillips l2003T Petfini & Pagel 2004). 
While these methods can be calibrated empirically at low 
metallicity, nebular photoionization models have provided the 
necessary calibration at high metallicity (around the solar 
value and above), because of the lack of reliable H II region 
abundances. 

From the observations of a few extragalac tic H II re gions 
that were believed to be extremely metal-rich (iKinkel & Rosal 
fT^lDTaTeTali riOOO: Castellanos, Diaz & Terlevich 2002; 
Kennicutt, Bresolin & Garnett 2003) a disconcerting discrep- 
ancy emerges between r P -based oxygen abundances and those 
derived through indirect means, the former being smaller by 
factors of a few. These results suggest that the oxygen abun- 
dance in the inner disks of spiral galaxies barely exceeds the 
solar value, while strong-line methods would suggest abun- 
dances 2-3 times higher. 

It is clear that this is exploration ground for large aperture 
telescopes, given the faintess of the auroral lines one needs 
to measure in order to obtain electron temperatures at high 
metallicity. A first step in this direction has been taken by Gar- 
nett, Kennicutt & Bresolin (2004), who used the 6.5m MMT 
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telescope on Mt. Hopkins, Arizona, to obtain electron temper- 
atures from the [Nil] A5755 line for two Hll regions in M51, 
CCM 10 and CCM 72. The Hll regions in this galaxy have 
long been considered to be extremely metal-rich, based on 
their very low excitation (weak [O III] A A4959,5007 lines rel- 
ative to H3) and on the results of photoionizatio n models by 
iDfaz et alJ (119911) . The metallicity obtained by iGarnett et al.l 
however, is roughly a factor of three smaller, barely 
exceeding the solar O/H value in the most metal-rich object, 
CCM 72. 

In this paper we present an enlarged sample of H II regions 
in M51, observed with the Keck I telescope. The deeper spec- 
tra we obtained allowed us to measure the [N II] A5755 line 
in ten H II regions, in conjunction with [S III] A6312 in seven 
of them. The larger number of objects and the measurement 
of two auroral lines increases our confidence about the con- 
clusions drawn in the previous work. These observations, to- 
gether with the data reduction, are described in § 2. The deter- 
mination of electron temperatures and chemical abundances is 
presented in § 3. Major results on the oxygen, nitrogen and 
sulfur abundances can be found in § 4. The impact of the re- 
sulting abundances on the calibration of strong-line methods 
is discussed in § 5. The Wolf-Rayet (WR) features detected 
among these H II regions are briefly presented in § 6. A few 
considerations on effective yields and the gas fractions in M5 1 
are addressed in § 7, and we summarize our conclusions in 
§8. 

2. OBSERVATIONS 

Spectroscopic observations of H II regions in M5 1 were car- 
ried out on the night of April 24, 2003 at the W. M. Keck Ob- 
servatory on Mauna Kea, with the Keck I telescope equipped 
with the Low Resolution Imaging Spectrometer ( O ke et all 
fl995l) . The night was photometric, with seeing at or below 
1". Using the blue and red channels of the spectrograph with 
a dichroic beam splitter, we obtained simultaneous observa- 
tions at wavelengths below and above ~ 5300 A. For the blue 
setting a 600 lines mm -1 grism blazed at 4000 A was used, 
providing spectra with a resolution of 4 A FWHM on a mosaic 
of two 2Kx4K Marconi CCDs down to the UV atmospheric 
cutoff. For the red spectra the main observations, covering 
roughly the 5300-7200 A wavelength range, were carried out 
with a 900 lines mm -1 grating blazed at 5500 A (3.5 A resolu- 
tion) and a 2Kx2K SITE CCD. Additional spectra with a 400 
lines mm -1 grating blazed at 8900 A (7 A resolution, 7000- 
10,000 A approximate coverage) were obtained for the mea- 
surement of the near-IR sulfur lines at AA9069,9532 A. How- 
ever, due to the limited near-IR extension of the calibration 
data, only the A9069 fluxes are useful, and the A9532 fluxes 
were estimated from the theoretical ratio A9532/A9069 = 2.44. 

Two slit masks were used for the multi-object spectroscopy, 
containing 21 and 16 l"2-wide slits, respectively, with slit 
lengths typically in the 10-20 arcsec range. The slits were 
oriented in the North-South direction, and the airmass during 
the observations was below 1.3. Archival V, R and Ha images 
of M5 1 provided the necessary astrometry for the mask prepa- 
ration. The target H II regions were selected among the bright- 
est available in M51. A finding chart is shown in Fig.[0 Five 
objects are in common with the previous work by Bresolin, 
Kennicutt & Garnett ( 1999): CCM 10, 53, 55, 71A. and 72 
(catalog numbers from Carranza. Crillon & Monnet 1969). 

The observations analyzed in this paper consist of 2 x 1 800 s 
exposures in each of the 600/4000 and 900/5500 setups for the 
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FIG. 1. — Identification of the Hll regions included in the Keck multi- 
object observations, plotted on top of a Kitt Peak 0.9m telescope Ha im- 
age of M5 1 (image courtesy T. Rector; N at top, E at leftV Numbers are 
from Carranz a. Crillon & Momiel H969I) . except for P 203 IPetit et all 19961) . 
Squares mark the H II regions analyzed in this paper. 



first mask (objects east of the galactic nucleus), and 3 x 1800 s 
in each setup for the second (west of the nucleus). Two 900 s 
exposures were used for the near-infrared (400/8900) spectra. 
The spectrum from each slit was treated as a separate spec- 
trum, and reduced with common IRAF 2 routines. Standard 
star spectra obtained at the start, middle and end of the night 
provided the flux calibration. A standard extinction curve for 
Mauna Kea was adopted for th e atmospheric extinction cor- 
rection (Rrisciunas et al. 1987). 

The spectral region in common between the red and near- 
infrared spectra was used for flux normalization, together 
with the flux in the Paschen 9 and Paschen 10 lines relative 
to H/3 when available. This procedure introduced some ad- 
ditional uncertainty in the [S III] A A9069,9532/[S III] A6312 
line flux ratio, from which the electron temperature in the 
intermediate-excitation zone is derived (see §3). Given the 
insufficient useful overlap between the blue and red spectra, 
emission line fluxes in the red spectra were adjusted relative 
to the fluxes in the blue spectra by scaling the flux in the Ha 
line so that Ha/H/3 = 3.00 as in case B at T e = 6,000 K. The ra- 
tio of the [N II] AA6548,6583 to [N n] A5755 line fluxes, from 
which the electron temperature in the low-excitation zone is 
derived, is unaffected by scaling errors, as these lines fall on 
the same spectra (the 900/5500 ones). Moreover, this ratio is 
also largely insensitive to uncertainties in reddening given the 
small wavelength baseline involved. 

The reddening coefficient C(H/3) was determined from the 
Balmer decrement using the H/3, H7 and US lines, adopting 
the case B theoretical rat ios at T e = 6,000 K, and the inter - 
stellar reddening law of ICardelli Clavton & Mathisl <1989tl . 
We have measured the [N II] A5755 auroral line in ten H II 
regions. Out of these objects, [S III] A63 12 has been mea- 
sured in seven cases. The IR spectrum of CCM 57A, how- 
ever, fell outside of the area covered by the CCD, and there- 
fore for this object we could not measure T[S III] in the way 
outlined in § 3. The [Olll] A4363 line remained undetected 

2 IRAF is distributed by the National Optical Astronomy Observatories, 
which are operated by the Association of Universities for Research in As- 
tronomy, Inc., under cooperative agreement with the National Science Foun- 
dation. 
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in all objects, while [Oll]A7325 was detected in CCM 10, 
CCM 57, CCM 71 A and CCM 72, however with a generally 
poor signal-to-noise ratio. Besides th e two nebulae in com- 
mon with the study by Garnettet al. (2004), this is the first 
time that auroral lines are measured in a large number of H II 
regions in M5 1 . We will concentrate the rest of our analysis 
on the ten objects with a measured [N II] A5755 line, as in- 
formation on electron temperatures and chemical abundances 
can be obtained directly only for these H II regions. The fluxes 
for their most important lines, corrected for reddening and 
normalized to F(H/3)= 100, are presented in Table 1, where 
each H II region is identified by its CCM number. The H II 
region P203, 7" West and 61" North of the M51 nucleus, not 
includ ed in the Carran za et al. (1969) catalog, is identified 
from iPetit et alJ d!996h . The errors quoted in the Table re- 
flect uncertainties in the flat fielding and the flux calibration, 
as well as statistical errors, which dominate in the case of the 
fainter lines (e.g. the auroral lines). The comparison con- 
cerning line fluxes for objec t s in co mm on with other a uthors, 
in particular Bresoli neTaTl i 19991) and Idarnett et all (120041) 
(CCM 10 and CCM 72), is generally good. The agreement 
with these works is typically at the 5 percent level or better for 
the [N II] and [S II] lines, and at the 15 percent level or better 
for the [Oil], [Oni] and [S III] lines. A notable exception is 
given by t he [O Hi] lines of CCM 53 (46 percent larger than 
in the lBresolin etalJl 999 paper), and the [S III] A A9069,9532 
lines of CCM 71A (a factor of two smaller). At least for 
the latter Hll region, a well extended object with multiple 
emission peaks, differences in pointing can be responsible 
for the disagreement. Large differences are also found for 
C CM 10 when com paring with the line intensities reported 
bv lDfazetalJ(IT99!h . 

3. ELECTRON TEMPERATURES AND ABUNDANCES 

As in our previous work on M101 (Kennicutt et al. 2003) 
and on the M51 Hll regions CCM 72 and CCM 10 
( iGarnett et alJ 12004). we adopt a three-zone model to de- 
scribe the ionization structure of an H II region. These zones 
are characterized by different electron temperatures T e , each 
referring to different, coexisting atomic ionization stages: 
[On], [Nil], [S ii] in the low-ionization zone, [Sin], [Arm] 
in the intermediate-ionization zone, and [O III], [Ne II I] in the 
high-io nization zone. The photoionization models of lGarnetH 
( 1992) predict simple scaling relations between the tempera- 
ture in the different zones, applicable in a wide T[0 III] range 
(2,000-18,000 K): 



T[S III] =0.83 T[Q III] + 1700 K, (1) 



T[NII] =T[OII] =0.70T[OIII] + 3000K. (2) 

Observational evidence in support of these relationships, in 
particular Eg. (1). , has been provided by a few authors (e.g. 
IGarnett et all 1 19971 IKennicutt etafl 120031) . One could con- 
sider alternative formulations, for exa mple the rela tionship 
between T[0 II] and T[0 III] given by llzotov et all ill 9941) 
based on the grid of photoionization models bv Ista sinska 
(1990). For the same measured value of T[Oll], their Eq. 
(4) would lea d to TTOin] lo wer by a few hundred degrees 
relative to the Garnett ( 1992) equation. However, if we limit 
the Stasinska models to T eff < 40,000 K for the ionizing stars, 



as seems appropriate for the H II regions in M51, we obtain 
good agreement betwen the two different sets of photoioniza- 
tion models. 

We established from the [S II] A6717/[S II] A6731 line 
ratio that all the Hll regions are in the low-density 
regime (N e < 150 cm -3 ). Subsequently, using the five- 
level atom program ne bular implemented in IRAF/STSDAS 
(|Sh aw & Dufounl 19951) . we determined the electron temper- 
atures T[N II] (ten objects) and T[S III] (six objects) from 
the emission line ratios [N II] AA6548,6583 / [N II] A5755 and 
[S in] AA9069,9532/[S m] A6312, respectiv ely. The results 
are su mmarized in Table 2. As in the M101 Kennic utt et alJ 
(2003) paper, we have updated the S III collisional strengths 
used b y nebular adopting the results by iTaval & Gupta 
(1999). 
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FIG. 2. — Relationship between T[NlI] and T[S III] m easured for Hll re- 
gions in M51 (this work: full dots) and M101 iKennicutt et al. 2003: open 
dots). The model prediction from Garnett 1 1992) is drawn as a dashed line. 

A good agreem ent is found with the independent T[Nll] 
determinations by IGarnett et all (|2004) for CCM 10 [7400 
(+1000,-600) K] and for CCM 72 (6000 ± 300 K). To verify 
the H II region zone model adopted for the following analy- 
sis, we have considered the correlation between T[N II] and 
T[S III] for the six M51 objects where both temperatures are 
available. As Fig.[2]shows, the observed relati on agr ees very 
well with the expectations based on the GarnettllfT9 92) models 
(dashed line). This result is encouraging, since it lends sup- 
port to the H II region zone model adopted here, indicating 
that the electron temperatures in the different excitation zones 
can be well established even when only one auroral line (e.g. 
[N II] A5755 at low excitation) is available. We have therefore 
calculated the 'missing' temperature for the high-ionization 
zone T[Olll], using both equations (1) and (2), and taking a 
weighted average for the result when both T[N II] and T[S III] 
are available. This gives us electron temperature estimates 
for the high-, intermediate- and low-ionization zones as de- 
scribed above; the results are listed in Table 3. For six of 
the ten objects listed we have two independent estimates of 
T[0 III], which we have averaged, from our measured T[N II] 
and T[S III] . We n ote that the 2a value of T[S III] » 5400 K 
for CCM 72 from Garnettet al. (2004) is very close to our 
result (5700 ± 300 K). The lowest temperature in our H II re- 
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TABLE 1 

Dereddened Line Fluxes and Errors 



Line CCM 10 CCM 53 CCM 54 CCM 55 CCM 57 CCM 57A CCM71A CCM 72 CCM 84A P203 



[Oil] 


3727 


126 ±9 


129 ±9 


115 ± 8 


78 ±5 


114 ±8 


104 ±7 


147 ± 10 


63 ±6 


125 ± 12 


32 ±2 


[Ne III] 3869 


0.7 ±0.1 


0.6 ±0.1 


1.1 ±0.2 


0.8 ±0.1 




12.8 ± 0.9 


2.2 ± 0.2 


0.4 ±0.1 


1.6 ±0.4 




H<5 


4101 


25 ±2 


24 ±2 


24 ±2 


23 ± 1 


24 ± 2 


24 ±2 


25 ± 2 


25 ± 2 


21 ±2 


21 ± 1 


Hi 


4340 


48 ± 3 


49 ±3 


46 ±3 


46 ±3 


47 ± 3 


48 ± 3 


46 ±3 


47 ± 3 


45 ± 3 


43 ± 3 


He I 


4471 


2.4 ± 0.2 


2.9 ± 0.2 


2.4 ±0.3 


2.1 ±0.1 


2.9 ± 0.2 


5.6 ± 0.4 


2.7 ± 0.2 


2.1 ±0.2 


4.8 ± 0.5 




[O III] 


4959 


4.4 ± 0.3 


11.8 ±0.7 


14.2 ± 0.9 


6.7 ± 0.4 


5.7 ± 0.3 


60.0 ± 3.5 


16.1 ± 1.0 


2.4 ±0.1 


29.8 ± 1.8 


1.7 ±0.2 


[0 III] 


5007 


12.1 ± 0.7 


33.3 ± 2.0 


41.5 ± 2.4 


19.0 ±1.1 


15.9 ± 0.9 


168.6 ± 9.9 


46.1 ± 2.7 


6.6 ± 0.4 


84.9 ± 5.1 


4.9 ± 0.3 


[Nil] 


5755 


.50 ±.04 


.54 ±.07 


.65 ±.07 


.43 ±.04 


.68 ±.09 


.48 ±.08 


.79 ±.07 


.28 ±.04 


.92 ±.15 


.15 ±.02 


He I 


5876 


7.7 ±0.5 


9.4 ± 0.6 


10.0 ±0.6 


8.4 ± 0.5 


10.1 ±0.6 


11.4 ±0.7 


9.3 ± 0.6 


5.6 ± 0.4 


13.1 ± 1.0 


4.5 ±0.3 


[Sin] 


6312 


.30 ±.03 


.31 ±.06 


.44 ±.06 


.31 ±.03 


.31 ±.07 


.62 ±.09 




.16 ±.03 






[Nil] 


6548 


37 ± 3 


40 ±3 


40 ±3 


38 ±3 


40 ±3 


23 ± 2 


39 ±3 


33 ± 3 


49 ±4 


26 ±2 


Ha 


6563 


300 ± 20 


300 ± 20 


300 ± 20 


300 ± 20 


300 ± 20 


300 ± 20 


300 ± 20 


300 ± 27 


300 ± 27 


300 ±20 


[Nil] 


6583 


112 ± 8 


120 ±8 


122 ± 8 


117 ± 8 


124 ±8 


70 ± 5 


116±8 


100 ±9 


144 ± 13 


79 ± 5 


He I 


6678 


2.0 ±0.1 


2.4 ± 0.2 


2.6 ± 0.2 


2.1 ±0.1 


2.3 ±0.2 


3.5 ± 0.3 


2.7 ± 0.2 


1.8 ±0.2 


3.1 ±0.3 


1.1 ±0.1 


[S II] 


6717 


27 ± 2 


24 ±2 


36 ±3 


23 ±2 


26 ±2 


20 ± 1 


41 ±3 


26 ±2 


19 ±2 


23 ± 2 


[S II] 


6731 


19 ± 1 


17 ± 1 


27 ±2 


18 ± 1 


19 ± 1 


14 ± 1 


28 ±2 


20 ±2 


15 ± 1 


18 ± 1 


[Arm] 7135 


2.2 ± 0.2 


2.7 ±0.3 


3.6 ±0.4 


2.6 ± 0.3 


2.6 ± 0.3 




3.9 ± 0.5 


1.4 ± 0.2 


8.8 ± 2.8 


1.7 ±0.2 


[S III] 


9069 


17 ± 3 


18 ±4 


17 ±3 


21 ±4 


14 ±3 




13 ±3 


14 ±3 


17 ± 3 


9±2 


C(H/3) 




0.25 ±0.05 


0.10 ±0.05 


0.15 ±0.05 


0.05 ±0.05 


0.15 ±0.05 


0.05 ±0.05 


0.25 ±0.05 


0.35 ±0.10 


0.60 ±0.10 


0.21 ±0.05 


EW(H/3) (A) 


51 


92 


144 


74 


155 


300 


96 


61 


203 


47 



TABLE 2 
Measured Electron 
Temperatures (K) 



ID 


T[NII] 


T[SIII] 




(K) 


(K) 


CCM 10 


6900 ± 200 


6400 ± 300 


CCM 53 


6900 ± 300 


6400 ± 400 


CCM 54 


7300 ± 300 


7200 ± 500 


CCM 55 


6600 ± 200 


6100 ± 300 


CCM 57 


7400 ± 300 


6800 ± 600 


CCM 57A 


7900 ± 400 




CCM71A 


7700 ± 200 




CCM 72 


6100 ± 200 


5700 ± 300 


CCM 84A 


7700 ± 400 




P203 


5600 ± 200 





TABLE 3 

Adopted Electron Temperatures (K) 



ID 


T(0+, N+, S + ) 


T(S +2 , Ar+ 2 ) 


T(0 +2 , Ne+ 2 ) 




(K) 


(K) 


(K) 


CCM 10 


6900 ± 200 


6400 ± 200 


5600 ± 200 


CCM 53 


6900 ± 300 


6400 ± 300 


5600 ± 300 


CCM 54 


7300 ± 300 


7200 ± 300 


6300 ± 300 


CCM 55 


6600 ± 200 


6100 ± 200 


5200 ± 200 


CCM 57 


7400 ± 300 


6800 ± 300 


6200 ± 400 


CCM 57A 


7900 ± 400 


7500 ± 400 


7000 ± 600 


CCM71A 


7700 ± 200 


7300 ± 200 


6800 ± 400 


CCM 72 


6100 ± 200 


5700 ± 200 


4600 ± 200 


CCM 84A 


7700 ± 400 


7300 ± 400 


6900 ± 600 


P203 


5600 ± 200 


4800 ± 200 


3700 ± 300 



gion sample is measured for P 203, close to the galactic center, 
which, at T[N II] = 5600 K, is the coolest extragalactic nebula 
to date where a direct electron temperature exists. 



The ionic and element abundances derived from the elec- 
tron temperatures in Table 3 using the five-level atom program 
are listed in Table 4 and 5, respectively. The associated errors 
derive mostly from the uncertainty in the electron tempera- 
tures. The uncertainty in the sulfur abundance is largest for 
CCM 57 A, where we had to rely on A6312 alone for the S +2 
ionic abundance, due to the lack of the 1R spectrum mentioned 
in § 2. 

Ionization correction factors for sulfur are expected to be 
small in low-ionization nebulae, as is the case here. This 
is verified in the top panel of Fig. |5J where the ion ratio 
(S + + S +2 )/(0 + + +2 ) is plotted against the oxygen ionization 
fraction + /0. The M5 1 data are in general agreement with 
the trend observed for Hll regions in M101 and NGC 2403 
(open circles and squares, respectively), extended to higher 
+ /0, even if a couple of objects (CCM 55 and 57) display 
a slight offset from the average relation. The dashed line 
in Fig. El reprod uces the s ulfur ionization correction formula 
(Stasinska 1 978tlFrenchll98l : 



where we have used a = 2.5 and log (S/(D)= -1.6, which are 
the parameters we adopted in our M101 paper. Our new ob- 
servations in M5 1 do not justify a change in the parameters 
that appear in the ionization correction formula. 

In the case of argon the situation appears more complex. 
Unlike sulfur, only one ionization stage (Ar +2 ) is observed 
for this atomic species in the optical spectra of H II regions, 
and in low-ionization nebulae the amount of unseen Ar + can 
become non-negligible. In our previous work on M101, we 
established that Ar +2 / S +2 is a good tracer of Ai7 S, based on 
the near independence of this ratio on the oxygen ionization 
fraction (at the same time we showed that Ar +2 / +2 is a bad 
tracer of Ar/O, and that neon has similar problems). The 
bottom plot of Fig. [3] suggests, however, that at high + /0 
(low ionization) Ar +2 / S deviates from the roughly constant 
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FIG. 3.— Top; The ion ratio (S + + S +2 ) /(0 + + +2 ) plotted against the 
oxygen fractional ionization + /0 for our sample of objects in M5 1 (full dots 
with error bars), and for objects studied in M101 i Kennicutt et al. 2003, open 
circles) and in NGC 2403 iGarnett et al. 1997, open squares). Equation (3) 
for the sulfur ionization correction factor provides the dashed curve. Bottom: 
Ar +2 /S +2 plotted against + /0 for the same sample of objects as above. A 
correlation with ionization seems to occur at large + /0 values. The vertical 
bars show the estimated correction factors for the M51 Hll regions (see text 
for explanation). 
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citation (as for most of the M101 and NGC 2403 H II regions 
plotted) a moderate negative correction (~ -0.1 dex) should 
be applied to the observed Ar +2 / S +2 ratios in order to recover 
the real Ar/S ratios. This will remove most of the appar- 
ent decreasing trend of Ar/S with oxygen ionization fraction 
suggested at first from Fig. [3] 

4. RESULTS FOR OXYGEN, SULFUR AND NITROGEN 

The somewhat surprising result we obtain concerns the 
typical oxygen abundance of the M5 1 H II regions stud- 
ied here: for most objects O/H is below the solar value 
[log (O/H) = -3.31, Allende Prieto et al. 2001], and reaches 
-1.4 times the solar value for P203 [log (O/H) = -3.16], 
the most oxygen-rich H II region in our sample. This con- 
trasts with the res ults of detailed photoionization models by 
iDiaz et all lll991ll . who estimated highly super-solar abun- 
dances, between log (O/H) = -2.9 and log (O/H) = -2.6 for 
their objects (CCM 10 and CCM 72 are in common with 
our sa mple), and thus confirms the findings of Garnett et al. 
(2004). The concept of M51 as a very metal-rich (largely 
over-solar) spiral galaxy, supported since the early 1990's, 
based on the above-mentioned models and the accepted cali- 
brations of several empirical abundance determination meth- 
ods, in particular R^i ( Vila-Costas & Edmunds 1992), is not 
fully supported by our data. The results for P 203, on the other 
hand, suggest that the central region of the galaxy, where the 
excitation is very low, contains nebulae with a O/H ratio at 
least 40% larger then the solar value. The discussion of ad- 
ditional low-excitation H II regions near the M5 1 nucleus is 
deferred to the end of § 5. 

Column 2 in Table 5 lists the deprojected distance from the 
galaxy center in units of the isophotal radius (Rq = 5.4 ar- 
cmin, corresponding to 13.2 kpc if the distance of 8.4 Mpc is 
adopted from Feldmeier et al. 1997). Our data span the disk 
from R /R = 0. 19 to R /R Q = 1 .04, but with the majority of the 
objects concentrated around R/Rq = 0.5. As a result the radial 
trends of abundances and abundance ratios depend heavily on 
the Hll regions at the two opposite extremes of the range, 
P203 (inner) and CCM 71A (outer). We plot the radial distri- 
bution of O/H, N/O and S/O in Fig.|5] which indicates a rather 
shallow O/H gradient in M5 1 . The dashed line represents the 
linear regression to the data points: 



FIG. 4. — The predicted correlation between the X(Ar + ')/X(S ) ionic frac- 
tion and the radiation softness parameter r\ from the models of Stasinska et 
al. (2001). The latter are plotted for an upper limit of the initial mass func- 
tion of 120 M (filled circles) and 30 M Q (open circles), for both extended 
and instantaneous bursts. The ionization correction factor used to obtain the 
abundance ratio Ar/S from Ar +2 /S +2 can be estimated from the reciprocal of 
X(Ar +2 )/X(S +2 ) in this plot. 



value found at higher ionization. This is likely an effect of 
the increase of Ar + with respect to Ar +2 at low excitation. 
There is, in fact, a good correlation between Ar +2 / S +2 and the 
hardness of the radiation field, as measured by th e parameter 
X] = (0 + /0 +2 )/(S + /S +2 ) dVflchez & Pagell 19881 see Fig.|4}. 
The ionization correction factor necessary to estimate Ar/S 
from Ar +2 / S +2 was obtained from the corresponding quanti - 
ties in the models of Stasinska, Schaerer & Leitherer (2001), 
and indicate that at low excitation the correction can amount 
to +0.2 to +0.3 dex. The shift for each H II region in the dia- 
gram thus determined is shown by the vertical bars in Fig. [3] 
This ionization correction scheme also shows that at high ex- 



12 + log(0/H) =8.72 (±0.09) -0.28 (±0.14) tf/flo (4) 

The gradient from equation (4) corresponds to -0.02 ±0.01 
dex kpc -1 . The central O/H abundance ratio suggested by the 
regression is about solar or slightly above solar, but we can- 
not exclude a steepening of the relation at small galactocen- 
tric radii. Clearly, data for additional H II regions are needed 
to better constrain the abundance gradient in M5 1 . There is 
considerable scatter at a given radial distance, amounting to 
—0.4 dex in O/H around R/Ro = 0.5. This is not atypical, 
since scatter of a similar magnitude is observed in other well- 
studie d spiral galaxies llKennicutt et al J 120031 IGarnett et alJ 
119971) . Finally, the radial trend of both N/O and S/O are con- 
sistent with a constant value, although a very shallow N/O 
gradient is suggested by Fig. [5] 

Plots of the abundance ratios S/O and N/O as a function of 
O/H are shown in Fig. [6] The S/O ratio (upper panel) for the 
M5 1 H II regions is consistent with the approximately con- 
stant value [log (S/O) =—1.6] observed in other spiral galax- 
ies, without any evidence for a decrease as the oxygen abun- 
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TABLE 4 
Ionic Abundances 



ID 


+ /H + 


+2 /H + 


N+/H+ 


S+/H+ 


S +2 /H + 


Ne +2 /H + 




Ar +2 /H + 


CCM 10 


(3.1±0.6)xl(r 4 


(5.5±1.4)xl(T 5 


(6.8±0.7)xl0~ 5 


(3.0±0.3)xl0~ 6 


(7.5±1.3)xl0~ 6 


(1.6±0.5)xl0 -5 


(7.2±1.5)xl(r 7 


CCM 53 


(3.1±0.8)xl(r 4 


(1.5±0.6)xl(r 4 


(7.3±1.2)xl0- 5 


(2.7±0.4)xl0 -6 


(7.9±1.9)xl0 -6 


(1.5±0.7)xl(T 5 


(8.9±2.6)xl(r 7 


CCM 54 


(2.1±0.5)xl(T 4 


(1.0±0.3)xl0- 4 


(6.2±0.9)xl0- 5 


(3.5±0.5)xl0- 6 


(5.6±l.l)xl(T 6 


(1.2±0.5)x 


10" 5 


(8.2±2.0)xl0~ 7 


CCM 55 


(2.6±0.6)xl(T 4 


(1.4±0.5)xl(r 4 


(8.5±1.2)xl0- 5 


(3.3±0.4)xl(T 6 


(l.l±0.2)xl0~ 5 


(3.2±1.4)x 


10~ 5 


(1.0±0.2)xl0" 6 


CCM 57 


(2.1±0.6)xl0 -4 


(4.4±1.8)xl0~ 5 


(6.1±1.0)xl(T 5 


(2.4±0.4)xl0~ 6 


(5.5±1.7)xl0~ 6 






(7.0±2.6)xl0~ 7 


CCM 57A 


(1.2±0.4)xl(r 4 


(2.6±1.4)xl(r 4 


(2.7±0.6)xlO- 5 


(1.5±0.3)xl0~ 6 


(6.2±3.0)xl0~ 6 


(8.2±5.6)x 


10~ 5 




CCM 71A 


(2.0±0.4)xl(r 4 


(8.1±2.4)xl(r 5 


(4.8±0.6)xl0- 5 


(3.3±0.3)xl0~ 6 


(4.1±0.5)xl0~ 6 


(1.7±0.6)x 


icr 5 


(8.4±l.l)xl(r 7 


CCM 72 


(4.1±1.3)xl(r 4 


(i.o±o.4)xi(r 4 


(9.6±1.5)xl0~ 5 


(4.8±0.7)xl0~ 6 


(9.0±1.9)xl0 -6 


(3.8±2.2)x 


10~ 5 


(7.2±1.8)xl(r 7 


CCM 84A 


(1.7±0.6)xl(r 4 


(i.6±i.o)xi(r 4 


(6.2±1.3)xl0~ 5 


(1.6±0.3)xl(T 6 


(5.3±l.l)xl0- 6 


(1.3±1.0)x 


10~ 5 


(1.9±0.5)xl(r 6 


P203 


(3.2±0.8)xl(r 4 


(3.6±2.9)xl(r 4 


(l.l±0.2)xl0^ 


(6.0±0.8)xl0~ 6 


(9.8±1.7)xl0~ 6 






(1.8±0.4)xl0~ 6 




— , -0.2 
^ -0.4 
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FIG. 5 .— Radial trends of 12 + log (O/H) (top), (b) log (N/O) (middle) and 
log (S/O) (bottom). The dashed line represents the linear fit to the O/H radial 
gradient (equation 4). 
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FIG. 6. — Top: The relationship between the S/O abundance ratio and O/H. 
Symbols are the same as in Fig. [3] A trend of approximately constant S/O is 
seen to extend beyond the solar oxygen abundance. Bottom: The relationship 
between N/O and O/H. Symbols as in Fig.l3l 



TABLE 5 
Total abundances 



ID 


R/R 


12+log (O/H) 


log (N/O) 


log (S/O) 


CCM 10 


0.59 


8.56±0.07 


-0.66±0.10 


-1.54±0.09 


CCM 53 


0.51 


8.66±0.09 


-0.63±0.14 


-1.63±0.13 


CCM 54 


0.56 


8.49±0.08 


-0.53±0.12 


-1.52±0.10 


CCM 55 


0.51 


8.60±O.O8 


-0.49±0.12 


-1.44±0.10 


CCM 57 


0.51 


8.40±0.09 


-0.54±0.14 


-1.50±0.14 


CCM 57A 


0.72 


8.58±0.14 


-0.65±0.17 


-1.61±0.24 


CCM71A 


1.04 


8.45±0.07 


-0.62±0.10 


-1.58±0.08 


CCM 72 


0.51 


8.71±0.10 


-0.63±0.15 


-1.56±0.13 


CCM 84A 


0.44 


8.52±0.13 


-0.44±0.19 


-1.64±0.17 


P203 


0.19 


8.84±0.16 


-0.46±0.13 


-1.61±0.20 



dance increases. The results of iDfaz et alJ flWfl) had pro- 
vide d some sug gestion that such a decrease might occur 
(e.g. lGarnettl2003l) . Based on our results, we therefore find no 
need to invoke changes in the massive star initial mass func- 
tion or in the nucleosynthesis at high metallicity (at least up 
to approximately the solar value). 



The nitrogen abundance, shown in relation to O/H in the 
bottom panel of Fig. [6] has an average value log (N/O) ~ 
-0.56, higher than the value measured in H II regions 
of comparable oxygen abundance in galaxies like M101 
and NGC 2403. This apparent nitrogen overabundance 
seems consistent with the known spread of N/O at a 
given O/H for spiral galaxies, related to differences in 
star formation histories, with earlier Hubble types (M51 
is classified as Sbc) h aving on average higher N/O ra- 
tios than later types ([Henry. Edmunds &' Kocoen 2000; 
iPilvugin. Thuan & Vilchezl2003l) . 

5. EMPIRICAL ABUNDANCE INDICATORS 

The results of a quarter century of research on nebular 
abundances have proved very useful to obtain estimates of 
abundance gradients in spiral galaxies and of oxygen abun- 
dances of star forming regions in the distant universe. Indi- 
rect methods (as opposed to the direct methods based on the 
determination of electron temperatures) using ratios of strong 
lines have been developed by several authors, the most widely 
used being arguably the R23 empirical abundance indicator of 
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iPagel et al.l lf!979). Some evidence has been mounting, how- 
ever, against the accuracy of such methods in the abundance 
range typical of the central regions of spiral galaxies (half so- 
lar oxygen abundance and above). This has been discussed at 
length in our work on M101 (Kennicutt et al. 2003), and we 
will not repeat the discussion here. In that paper we also pre- 
sented the main issues which could complicate the analysis 
based on the collisionally excited lines, namely the presence 
of temperature fluctuations and, at high metallicity, temper- 
ature gradients. Their effect is similar, in that the measured 
electron temperature is weighted towards the higher temper- 
ature zones, leading to a possible underesti mate of the real 
nebula r abundances. The reader is referred to Kennic utt et al.l 
(2003) for further details and for a justification for the use of 
the collisionally excited lines. Admittedly, the issue is still 
open. In addition, the temperature in the [Nil] zone could 
be overestimated from neglec ting the recomb ination compo- 
nent of the A5755 line (see iLiu et alJl2000h . The lack of 
the line flux from [N III] at 57 for our M5 1 sample does 
not allow us to estimate the contribution of this recombina- 
tion directly from observations. However, this contribution 
is likely to be small, given the small amount of N" 1 " 1 " present 
in low-excitation nebulae. The good agreement between the 
obse rved T[Nll]-Tr S III] relationship and the predicted one 
from lGarnettl (11992) also suggests that such a recombination 
component can be neglected to first order. 

The detection of auroral lines in objects around the solar 
O/H ratio is prompting a revision of the abundances based 
on indirect methods, since relatively low abundances are 
found for objects that were once th ought to be very metal- 
rich (i.e. well above the solar value: Castel lanos et alJl2002l 
iKennicutt et alJl2003l) . This is also the case for the M51 H n 
regions analyzed in this work, and a brief assessment of a 
few representative empirical abundance indicators seems to 
be warranted at this point. 

It is obvious that the abundances obtained in M51 in the 
present work have an important impact on most empirical 
calibrations. This is exemplified in Fig. Q where we show 
the comparison of our direct oxygen abundances with those 
obtained by means of j?23, adopting two dif ferent calibra- 
tio ns: the one by [Edmunds & Pa gefl (119841) and the one 
by iPilvuginl ( 1200 ll) (see iGarnett et al.l 120041 for a compari- 
son involving two other calibrations). The sample of ob- 
jects shown includes, besides the M51 Hll regions stud- 
ied here, a number of nebulae extracted from the litera- 
ture (Kennicutt et al. 2003; Garnett et alJI 19971 Ivan Zee et al.l 
1 19981 lDeharvengetalJl2000t iDfaz & Perez-Monterol l2000t 
IDfaz et alJl2000t ICastellanos et al.ll2002l) . We note that for 
the high-metall icity object CDT1 in NGC 1232 studied by 
ICastellanos et alJ d20 02) we derive, using their line fluxes, an 
O/H value smaller by 0. 1 3 dex than their published abundance 
12 + log (O/H) = 8.95. Our results on the M51 Hll regions 
clearly demonstrate the failure of the existing R23 calibrations 
at low excitation (log7?23 < 0.5), since they overestimate the 
abundances by 50^-00 percent. This discrepancy corresponds 
to electron temperatures smaller by approximately 10-15 per- 
cent compared with our measured values (about 3a to 4a). In 
the particular case of CCM 72 and CC M 10, the [Oil] tem- 
peratures modeled bv lDfaz et al] (119911) (4800 K and 6000 K, 
respectively) lie 4a and 6a below our measurements. 

This result has fundamental consequences for the studies of 
abundance gradi ents in spiral galaxies based on this em pirical 
indica tor (e.g. IZaritskv et al.lll994l IVila-Costas & Edmunds 
1992), and the actual gradients in the central regions will in 
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FIG. 7. — Comparison between O/H abundances obtained from measured 
electron temperatures (dots) and from the semi-empirical abundance indi- 
cator R23 (full and dotted lines). Two different R23 calibrations are used: 
Edmunds & Pagel (1984, = EP, full line) and Pilyugin (2001, = P01, dotted 
lines). Two curves are plotted for the latter, corresponding to two different 
values of the excitation parameter, P = 0. 1 and P = 0.9. The plotted data are 
from t h e current work (full dots) and from published stu dies: Kennic utt et a ll 
l2003l).lGarnett et alJ Il99l|) I van Zee et alJ <199l3).|Deharyeng et alJ 1200(1) 
Diaz & Perez-Montero 1 2000) and Castellanos et al. 1 2002) (open squares). 
The lowest excitation object, a t the far left, is P203, w hile the most metal- 
rich one is CDT1 in NGC 1232 1 Castellanos et al. 2002). The dashed line is a 
fit to the M5 1 points, except CCM 84A: 12 +log (O/H) = 8.64-0.33 logR 23 . 



general be much flatter than reported in these works. The 
discrepancy is also n ot peculiar to /?23- As an example, 
Pettini & Pagel (2004) have considered an empirical abun- 
dance calibration of two line indi ces involving the ni trogen 
lines: [Nil] A6583/Ha (see also iDenicolo et alJl2002T) . and 
([Oiii] A5007/H/3)/([N 11] A6583/Ha). We show in Fig.[8]the 
position occupied in the corresponding diagrams by the M5 1 
Hll regions. A decrease in the slope at low exc itation is no- 
ticed with respect to the Pettini & Pagel (2004) calibrations 
(shown by the dashed lines), indicating a flattening around 
12+log (O/H) = 8.6, which however also reflects the large ni- 
trogen abundance found in M51 from our analysis. Fig. [5] 
shows inst ead how the O/ H abundances derived through the 
P-method (Pilvugin 2001) compare with direct abundances. 
In this plot, only objects which are predicted by the P-method 
to have 12 + log(0/H)>8.2 (the quoted range of validity for 
this method in the high-abundance regime) are plotted. This 
indicator seems to be afflicted by severe difficulties in a wide 
abundance range. The general conclusion we can draw is that 
great caution is necessary when using strong-line empirical 
abundance indicators for low-excitation H II regions. Erro- 
neous results can be obtained in the inner disks of spirals up 
to a factor of a few. 

How low can the excitation in M51 H II regions be? We re- 
port in Table 6 the line fluxes, normalized to H/3 — 100, mea- 
sured for the H II regions belonging to the central zone of the 
galaxy (R = 0.12-0.36 R ) besides P203. No auroral line 
was detected in these objects, therefore not permitting a di- 
rect abundance determination. The total [Olll] AA4959,5007 
flux is just a few percent of the H[3 line flux, and the corre- 
sponding value of log/?23 is m the range -0.42 to -0.67. Be- 
cause of the rather flat relationship between O/H and ^23 at 
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FIG. 8. — Top: O/H abundances obtained with the direct method are plotted 
against t he line ratio [Ni l] A6583/H a. Data are from thi s wo rk (full dots) 
and fr om Kennicutt et al. 1 2003 ), Garnett et al. 1 1 997 ) and Castella nos et alj 
12002]) (open circles) . The dashed line is the regression taken from 
Pettini & Page] 1 2004), obtained from a fit to a larger sample of H II regions. 
Bottom: Same as above, for the ([O III] A5007/H/3)/([N II] A6583/Ha) index. 




We have estimated upper limits for the [N II] electron tem- 
perature and the corresponding lower limits for the oxygen 
abundance from the non-detection of the [N II] A5755 line in 
the spectra of the central objects. The results, summarized in 
Table 6, are not very compelling, and they are still consistent 
with abundances in the central region of M5 1 around the value 
found for P 203. 

We finally point out that the He I A5876/H/3 line ratio cor- 
relates with both O/H (for the objects in Table 1) and galac- 
tocentric distance (objects in Tables 1 and 6), as shown in 
Fig- El Together with the reduced excitation measured from 
the [O lll]/H/3 ratio mentioned above, Fig. ^| shows how 
the ionizing field becomes softer towards higher metallicities, 
above a certain abundance threshold (the HelA5876/H/3 ra- 
tio decreases below its saturation level, at which He is com- 
pletely singly ionized, for 12 + log (O/H) > 8.5, top panel) and 
towards the galactic central regions (lower panel). The ob- 
served increase of the neutral helium fraction towards the 
metal-rich central region of M5 1 can be viewed as the result 
of the decrease in the equivalent effective temperature of the 
ionizing radiati on field wit h increasing metallicity (see also 
iBresolin & KennicuttfeOOH) . 
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FIG. 9. — Comparison betw een direct O /H abundances and those ob- 
tained through the P-method of Pilvugin 1 200 11). The data are from t he cur - 
rent work (full dots) and from Kennicutt et al. 1 2003), Garnett et al. 1 1997), 
Ivan Zee et all 1 1998), Deharveng et al. 1 2000), Diaz & Perez-Montero] <200Cft 
and Castellanos et al. 1 2002 ) (open squares). Only objects which are pre- 
dicted by the P-method to have 12 + log (0/H)>8.2 are plotted. 



low excitation suggested by the M51 points in Fig. [7] these 
very low R23 values might still correspond to metallicities 
not much higher than solar [12 + log (O/H) w 8.8-8.9], how- 
ever any extrapolation to determine an oxygen abundance is 
at present still very uncertain. The dashed line in Fig.0shows 
such an extrapolation, determined from the M51 H II regions. 
The highest oxygen abundance in our sample, correspond- 
ing to log/? 2 3 = -0.66 (the value for CCM 6A, CCM 37A 
and CCM 107 within ±0.01), would be 12 + log (O/H) w 8.86 
(1.5 times the solar O/H), comparable to our measured oxygen 
abundance for P 203. 



FIG. 10. — Top: A trend of decreasing He I A5876/H/3 line ratio with abun- 
dance above 12 + log (O/H) f» 8.5 is suggested here. The A5876 line has been 
corrected for an underlying absorption equivalent width of 0.8 A. Objects 
plotted are from M51 (filled dots) and M101 (open circles). Bottom: Varia- 
tion of He I A5876/H/3 with galactocentric distance for M5 1 Hll regions. The 
filled dots are objects with a direct abundance determination, the open stars 
are objects of low excitation found in the central region of the galaxy. 



6. WOLF-RAYET STARS 

The main spectroscopic signature of WR stars (the blue 
bump at 4660 A) has been detected in six of the H II regions 
examined in this work. As shown in Fig.^2 we can clearly 
distinguish the stellar NmA4634-41 and HellA4686 lines 
in all cases, together with CIVA4658. Given the absence of 
N V A4603-20, the morphology of the blue bump is character- 
istic of late WR subtypes (WNL). In addition, the C III A5696 
line from WC subtype stars has been clearly detected in the 
two most metal -rich Hll regions, CCM 72 and P203. WR 
features have been detected in other M5 1 H II regions not an- 
alyzed in this paper. In CCM 13, 24, 41, 84 A and 91 only 
the blue bump is present in our spectra. Both the 4660 A and 
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TABLE 6 

Line fluxes of low-excitation objects 



ID 


R/Ro 


C(H/3) 


EW(H/3) 


[On] 


[OlII] 


He I 


[Nil] 


[SII] 


T[Nn] 


12 + log 


(O/H) 








(A) 


3727 


4959+5007 


5876 


6548+6584 


6717+6731 


(upper limit) 


(lower 


limit) 


CCM6A 


0.12 


0.10 


55 


17±1 


4.4±0.3 


3.9±0.3 


80±4 


26±1 


5800 K 


S 


.40 


CCM 37A 


0.16 


0.09 


12 


17±1 


4.8±0.4 


3.0±0.3 


54±3 


20±1 


6200 K 


8 


.20 
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the 5696 A features are seen in CCM 37 A, 84 and 107, while 
only the 5696 A line is visible in CCM 6 A and 71 (also with 
5810 A). A marginal detection of He IIA4686 has been made 
in a few additional objects, as well. Note that all the nebulae 
(except CCM 7 1 ) where the C III A5696 line has been detected 
are located at small galactocentric radii, or have a large O/H 
abundance (CCM 72). This agrees with the known trend of in- 
creasing fracti on of WC relative to WN stars with abundance 
(Mass evl2003l) . 

The number of WNL stars can be estimated from the to- 
tal He II A4686 line luminosity, assuming a fixed stellar lu- 
minosity of 1.6 x 10 36 ergs" 1 dSchaerer & Vaccalfl998T) . At 
the adopted distance to M51 of 8.4 Mpc the two objects with 
the highest He II luminosity (CCM 10 and CCM 72) are then 
found to contain five-six WNL stars each. The luminosity of 
the remaining objects is consistent with just one or two WNL 
stars being present. These numbers are not particularly re- 
markable when compared, for exampl e, with the very large 
population of WR stars discovered by Crowther et al. (2004) 
in H II regions of another supposedly metal-rich galaxy, M83. 
Reporting the detection of WR features in solar-like abun- 
dance environments, however, is still of considerable inter- 
est, because of the trends one can establish between WR 
star population properties (e.g. their subtype distribution) and 
metallicity. It is also relevant to point out that WR fea- 
tures are found in very cool nebulae, such as the metal-rich 
CCM 72 and P203, in support of the idea that the pres- 
ence of these evolved massive stars does not significantly 
affect the ionizing propertie s of the embedded star clusters 
jBresolin & Kennicutl2002h . 

7. GAS FRACTIONS, EFFECTIVE YIELDS, AND CHEMICAL 
EVOLUTION 

A comparison of observed abundances and gas fractions in 
galaxies provides interesting information on the effect of gas 
flows on galaxy evolution. Variations of the effective yield, 
defined as 



where Z is the metallicity and /i is the g as mass fraction, 
are sensitive to gas inflo w and outflow (Edmunds 119901 
Koppen & Edmunds II 19991) . iGarnettl J2f)02) compared effec- 
tive yields, derived from global abundance, stellar mass, and 
gas mass properties, in a sample of 50 nearby spiral and ir- 
regular galaxies, and found that y e g varies systematically with 
galaxy rotation speed (and thus presumably total mass). The 
adopted interpretation was the low-mass irregular galaxies, 
with smaller effective yields than those derived for massive 




4600 4700 4800 

Wavelength (A) 



FIG. 1 1 . — The WR stellar lines around 4660 A have been detected in these 
six Hll regions out of the ten analyzed. Feature identification is provided by 
the vertical dotted lines. 



spirals, were losing a large fraction of their metals to the in- 
tergalactic medium (IGM). 

Based on our new abunda nces fo r M51, we can re-evaluate 
its effective yield. Garnett (2002) used the measured abun- 
dance at one effective disk radius as a measure of the average 
metallicity of a spiral disk. From his Table 4, the adopted av- 
erage 12 + log (O/H) was 9.12 for M51. Our new value for 
12 + log (O/H) at /4ff = 0.6 fl is 8.55, a factor of 3.7 smaller. 
Thus the new value of y e ff w e would derive is smaller by the 
same amount. TGarne tt (2002) derived y e ff = 0.012 for M51, so 
we would obtain y e g = 0.0032 from our new data. This lower 
value is comparable to the values derived for smaller irregular 
galaxies such as IC 10 and NGC 6822 in Garnett's sample. 
Note also that this value of y e g is much closer to that expected 
from stellar nucleosynthesis, and to the value obtained from 
estimates of the mean metallicity and gas fraction in the so- 
lar neighborhood. If other spirals show similar reductions in 
y e ff when more temperature-based abundance measurements 
become availabl e, it may b e necessary to re-examine the con- 
clusions of Garnett ( 2002) regarding what galaxies may expe- 
rience significant loss of metals. 

It is also interesting and important to look at the variation 
of y e ff across a disk galaxy, to see if gas flows may play a 
role in the chemical evolution within the disk. A comparison 
between gas fraction and abundance across the disk has been 
done for only a few spirals to date. 
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We obtained the radial distribution of gas surface densi- 
ties in M5 1 from the single-dish CO mapping of iKuno et al.l 
i ll 9951) . m easured with a 16" beam, and from the H i 21- 
cm map of iTilanus & Alien! (119911) . lYoung et all i ll 9951) also 
measured the radial distribution of CO, but their much larger 
beam size means that the mass surface density is averaged 
over a much larger area. We c hoose not to use m easurements 
based on interferometers (e.g. lAalto et al.lll9 99), as they can 
miss significant amounts of extended emission due to miss- 
ing short an tenna baselines; we note that the measurements of 
lAalto et all (119991) show a very s harply peaked CO flux distri- 
bution in rough agreement with Kuno et al. ( 1995). The de- 
rived molecular gas surface densities depend critically on the 
conversion from CO intensity to H2 column density. We de- 
rived molecular surface densities under two assumptions for 
the I(CO)-N(H?) conversion f actor: (1) the 'standard' factor 
3.0 x 10 20 cm" 2 (Kkm s" 1 )" 1 (Wils on! 19951) . which includes a 
contribution for heliu m; (2) the smaller f actor 1 .3 x 10 20 cm -2 
(K km s -1 ) -1 used bv lKuno et alJ (11995) (adjusted here to in- 
clude helium). Whichever value is used, the mass of molecu- 
lar gas in M5 1 is quite large, and the mass surface density of 
H2 exceeds that of atomic gas over much of the disk. The CO 
measurements do not go beyond 180" from the nucleus, so we 
have smoothly extrapolated the H2 surface densities to larger 
radii based on an approximate exponential fit to the inner disk 
profile. Gas surface density distributions for a subset of our 
assumptions are shown in Figure 12. 




R(arcsec) 

FIG. 1 2. — Gas and stellar mass surface density distributions in M51. The 
dashed line shows the H I surface densi t y (a djusted for the presence of he- 
lium) from the data of Tilanus & Allejj 119911) . The dotted lin e shows the 
molecular gas profile obtained from the CO measurements of K uno et aT] 
1 1995), assuming a CO-H2 conversion factor of 3 x 10 20 cm -2 (K km s -1 ) -1 . 
The dot-dash line shows the mass profile for the stellar exponential disk deter- 
mined as discussed in the text. The filled circles show stellar mass densities 
derived directly from the surface brightness data of Kuchinski et al. 1 2000). 

The mass surface density for st ars in M5 1 was derived from 
the BVRI surface photometry of lKuch inski et alJ ( 120001) . We 
converted the I-band surface brightnesses to mass surface den- 
sities by estimating the mass-to-light ratio (M/L) from the re- 
lations between M/L and c olor de rived from stellar population 
models in Bell & de Jong (2001). From the Kuchinski et al. 
surface brightnesses, we estimated a central I-band disk sur- 
face brightness of 19.05 mag arcsec -2 and a disk scale length 
of 1 14"; we adopted a linear color profile with radius rang- 
ing from B-R = 1.2 at the center to B-R = 0.9 at 300". The 
resulting stellar mass profile for the disk in M5 1 is shown in 
Figure 12 as the dot-dash line, while the mass densities es- 
timated directly from the photometry are shown as the filled 




100 200 300 

R(arcsec) 



FIG. 13. — (a) The radial variation of the gas mass fraction for three dif- 
ferent assumptions. Solid line: molecular gas masses estimated using I(CO)- 
N(H2) conversion of 3. Ox 10 20 ; dashed line: molecular gas mass densities 
estimated using I(CO)-N(H2 ) conversion of \3xi0 20 ; dotted line: same gas 
densities as for solid line, but with stellar surface brightnesses adjusted for 
extinction A(I) = 0.5 mag. (b) Radial variation of the gas fraction, where 
the stellar mass densities were derived directly from the I-band photometry, 
rather than the underlying exponential disk. Solid line: molecular gas den- 
sities obtained using I(CO)-N(H2) = 3.0xl0 20 ; dotted line: molecular gas 
densities obtained using I(CO)-N(H 2 ) = 1.3x 10 20 . 



circles. Note that the photometry shows evidence for a central 
spheroidal component. 

Inspection of Figure 12 shows one curious feature: the gas 
surface density profile is steeper than the stellar disk profile. 
This implies that the gas fraction decreases radially outward. 
We illustrate this in Figure 13(a), which shows the radial vari- 
ation in three different cases: (1) I(CO)-N(H 2 ) = 3.0 xlO 20 
(solid line); (2) I(CO)-N(H 2 ) = 1.3xl0 20 (dashed line); and 
(3) a case similar to (1), but in which the stellar surface bright- 
nesses have been adjusted upward for an extinction A(I) = 0.5 
mag (dotted curve). In all three cases, the derived gas fraction 
declines by a factor 3-4 radially outward. This behavior con- 
trasts strongly with that seen in other spiral gal axies, where 
the ga s fraction increases outw ard (M 81 - iGarnett & Shield's! 
I987t NGC 4254 -iHenrv et"allfl994l NGC 2403 and M33 - 
Garnett et all 19971) . 

We have also looked at what happens if we simply use the 
observed surface brightness profile to derive the stellar mass 
densities, rather than an underlying exponential disk. The re- 
sulting gas profiles are shown in Figure 13(b), for our two 
assumptions of the I(CO)-N(H2) conversion. While the re- 
sulting gas fractions for the inner disk are certainly smaller, 
due to the high central surface brightness, we still see a de- 
clining gas fraction with radius for R > 50". 

We note that a radial variation of a factor three in the CO- 
H 2 conversion is not sufficient to mitigate the decline in gas 
fraction. A factor six change in the conversion factor would 
be needed to flatten the gas fraction profile; such a variation 
is not at all predicted even by the most extreme estimates of 
metallicity variation in the conversion factor i llsraell997l) . A 
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radial decline of about 1 .5 mag in A(I) from the center to the 
outer d isk would also be en ough to flatten the gas fraction 
profile. IScovilleet all J200H) found a rough radial variation in 
A(V) for H II regions in the disk of M51, from A(V) « 4 in 
the inner disk to A(V) « 2 at R = 300"; however, it is not clear 
how to relate the extinction for H II regions, which are asso- 
ciated with gas and dust, to the extinction for the stellar disk. 
We also note that A(I) = 0.5 corresponds to E(B-V) = 0.3; 
with an observed B-V of 0.6, such a correction for reddening 
would make M5 1 a very blue galaxy unless the reddening is 
very gray. 

A radially decreasing gas fraction combined with a radi- 
ally decreasing metallicity implies an effective yield that in- 
creases with metallicity. Such a variation was inferred by 
IVila-Costas & Edmunds! dl992l) . but with older abundance 
data derived mainly from strong emission line calibrations. 
On the other hand, galaxies with abundances based on elec- 
tron tempe ratures show little ev idence for variations in effec- 
tive yields JGarnett et alJl997l) . 

The peculiar gas fraction profile in M5 1 may also point to 
dynamical influences. M51 is an interacting system. Fly- 
by interactions of this type can lead to bar formation, which 
induces radial flows of gas toward the center of the galaxy . 
iPiercd (QUI, iKohno et alJ II996T) . and lAalto et all 1i993) 
have presented photometric and kinematic evidence for a bar 
in the central regions of NGC 5194, while infrared imaging 
reveals that NGC 5195 has an obvious bar (e.g., the I-band 
images of Kuchinski et al. 2000). It is possible that the inter- 
action and bar have driven gas into the central regions of M5 1 , 
leading to the very steep gas profile we see. This would com- 
plicate the interpretation of the chemical evolution of this sys- 
tem, and requires dynamical and photometric modeling that is 
beyond the scope of this work. 

8. CONCLUSIONS 

We have analyzed a sample of ten Hll regions in the 
spiral galaxy M51 where at least one of the auroral lines 
[Nil] A5755 and [S III] A6312 could be measured, and we dis- 
cussed the resulting oxygen, sulfur and nitrogen abundances. 
The single, most remarkable result obtained is that in M5 1 the 
direct O/H abundances are considerably bel ow the value ob- 
tained either from photoionization models (iDfaz et al.ll 19911) 
or from empirical abundance indicators. The calibration of 
the latter depended so far, in the low-excitation range defined 
mostly by Hll regions in M51, on the abundances obtained 



from the models, therefore the two are not independent. We 
have not attempted to generate improved nebular models, for 
example including more recent stellar atmospheres for the 
treatment of stellar ionizing fluxes. This aspect needs to be 
explored in the future. 

The result concerning O/H has important implications 
for the calibration of empirical abundance indicators, and 
in general for the determination of abundance gradients 
in spiral galaxies. Our result, combined with previous 
ones obta ined for smaller sample s of metal-r i ch Hl l re- 
gions bv ICastellanos et alJ J2002I) . IDfaz et alJ fcOOOl) and 
Kennicutt et al. (2003), points to shallower gradients. Hll 
regions once believed to have O/H abundances equal to 2-3 
times solar are found by the direct method to rarely exceed the 
solar value, and if so by a moderate amount, up to 50 percent. 

Our analysis has also established that in our M51 sample 
the abundance ratio S/O is similar to the value found at lower 
O/H abundance in other spiral galaxies. We find no evidence 
for a decrease in S/O with increasing O/H, which would occur 
in case of a differing initial mass function and/or nucleosyn- 
thesis for massive stars. The nitrogen abundance has a mean 
value log (N/O) ~ -0.6, larger than in later spirals like M101, 
but this can still be interpreted with the known spread of N/O 
at a given oxygen abundance. 

The lower oxygen abundance we measure in M5 1 with re- 
spect to previous determinations has also allowed us to revise 
the effective yield of this galaxy. From a rather peculiar gas 
fraction radial profile we infer an effective yield increasing 
with metallicity. 

The number of galaxies with reliable direct measurements 
of nebular chemical abundances extending over their whole 
disks, including the metal-rich central zones, is still very lim- 
ited. These measurements are necessary for the correct de- 
termination of the chemical compositions in galactic disks 
and their interpretation in terms of chemical evolution mod- 
els. The current work on M5 1 demonstrates that large samples 
of metal-rich H II regions can be efficiently studied in spiral 
galaxies with current telescopes and instrumentation. 
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